NAS A Technical Memorandum 106629 
AIAA-94-2079 


Experimental Investigation of Jet-Induced 
Mixing of a Large Liquid Hydrogen 
Storage Tank 


C. S. Lin 

Analex Corporation 
Brook Park, Ohio 

M. M. Hasan and N. T. Van Dresar 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

6th Joint Thermophysics and Heat Transfer Conference 
cosponsored by the American Institute of Aeronautics and Astronautics 
and the American Society of Mechanical Engineers 
Colorado Springs, Colorado, June 20-23, 1994 



National Aeronautics and 
Space Administration 


( NASA-TM-106629) EXPERIMENTAL 
INVESTIGATION OE JET-INDUCED MIXING 
OF A LARGE LIQUID HYDROGEN STORAGE 
TANK (NASA. Lewis Research Center) 

12 D 


N9A-33607 


Unci as 


G3/34 0012113 


EXPERIMENTAL INVESTIGATION OF JET-INDUCED MIXING OF A 
LARGE LIQUID HYDROGEN STORAGE TANK 


C.S. Lin 

Analex Corporation 
3001 Aerospace Parkway 
Brook Park, Ohio 44142 

M.M. Hasan and N.T. Van Dresar 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Abstract 


Ja Jakob number, Ja = C p (T s - T b ) / 


Experiments have been conducted to investigate the 
effect of fluid mixing on the depressurization of a large liquid 
hydrogen storage tank. The test tank is approximately ellipsoi- 
dal, having a volume of 4.89 m 3 and an average wall heat flux 
of 4.2 W/m 2 due to external heat input. A mixer unit was 
installed near the bottom of the tank to generate an upward 
directed axial jet flow normal to the liquid-vapor interface. 
Mixing tests were initiated after achieving thermally stratified 
conditions in the tank either by the introduction of hydrogen 
gas into the tank or by self-pressurization due to ambient heat 
leak through the tank wall. The subcooled liquid jet directed 
towards the liquid- vapor interface by the mixer induced vapor 
condensation and caused a reduction in tank pressure. Tests 
were conducted at two jet submergence depths for jet Reynolds 
numbers from 80 000 to 495 000 and Richardson numbers 
from 0.014 to 0.52. Results show that the rate of tank pressure 
change is controlled by the competing effects of subcooled jet 
flow and the free convection boundary layer flow due to 
external tank wall heating. It is shown that existing correla- 
tions for mixing time and vapor condensation rate based on 
small scale tanks may not be applicable to large scale liquid 
hydrogen systems. 


Nomenclature 


A s liquid-vapor interface area 


k thermal conductivity 

m c average mass condensation rate at the liquid-vapor 

interface 


P, Pj, P f tank pressure, initial tank pressure, and final tank 
pressure 

Pr Prandtl number of bulk liquid, v/a 


QbiW 


volume flow rate of the free convection boundary 
layer along the tank wall 


Qc(x) 


jet volume flow at a distance of x downstream from 
the jet nozzle 


Qb.* 

Qc* 

Q j 

% 

Ri 


Q bl (x) evaluated at x = H| 

(^(x) evaluated at x = Hj 

jet volume flow rate at the jet nozzle outlet 

wall heat flux 

Richardson number evaluated at the begining 


each test, gP(T s -Tj) i Hj / Uj 2 


of 


C p specific heat of liquid 


Rej jet Reynolds number, Uj d/v 


D, D s tank diameter, liquid- vapor interface diameter 

d jet nozzle diameter 

g gravitational acceleration 


St c average condensation Stanton number, St c = m c 

V [A sP C p V s( T s- T b)l 

T b , Tj, T s bulk liquid temperature, jet temperature, and 
interace temperature 


Gr* the modified Grashof number, Gr* = g(Jq w H 4 / k v 2 T dome dome temperature 


Hj, Hj liquid height, jet submergence depth (T s -T b ) liquid subcooling 

h fg latent heat of condensation (T s -Tj) jet subcooling 
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V l m* 

mixing time, dimensionless mixing time 

U j 

jet velocity 

V . 

liquid volume 

v s 

average r.m.s turbulent velocity (or turbulence 
intensity) at the liquid-vapor interface 

X 

height (vertical distance) measured from the tank 
bottom or downstream distance from 
the jet nozzle 

Greek Symbols 

a 

thermal diffusivity of liquid 

p 

thermal expansion coefficient of liquid 

V 

kinematic viscosity of liquid 

P 

liquid density 

a 

surface tension 


Introduction 

Pressure control of a cryogenic storage tank in low 
gravity is one of the key advanced technologies being devel- 
oped in the cryogenic fluid management program at NASA. In 
pressure control technology, fluid mixing such as axial jet- 
induced mixing is considered an efficient method to reduce the 
pressure of a cryogenic tank. In normal gravity, when a 
partially filled cryogenic tank is subjected to external heating 
and the fluid is thermally stratified, the tank pressure increases 
with a rate usually higher than its homogeneous pressure rise 
rate 1 . Due to the temperature stratification the bulk liquid is 
generally subcooled. By using a mixer located near the bottom 
of the tank, the tank liquid is circulated and the subcooled axial 
jet can be directed towards the liquid- vapor interface. The 
subcooled jet spreads from the jet nozzle outlet, but remains 
confined to the axial region. The jet impinges on the liquid- 
vapor interface, which is kept horizontal by gravity, turns 
radially outward, and eventually deflects down into the bulk 
liquid near the tank wall and mixes with the bulk liquid. This 
mixing of the tank liquid by the jet destroys the liquid tempera- 
ture stratification and induces vapor condensation at the 
interface, resulting in the reduction of tank pressure. The 
required mixing time and vapor condensation rate at the 
interface are the two key elements in the mixing process to 
reduce tank pressure. The high jet flow rate fluid mixing and 
the associated vapor condensation results, which are not 
influenced by the buoyancy effects, may be applicable for the 
design of space-based systems. 

Use of the jet-induced mixing process to control the 
thermodynamic state of propellents in storage tanks has been 


extensively investigated. Several mixing time correlations 2 " 6 
have been proposed based on the concentration equilibrium, 
thermal equilibrium, or pressure equilibrium. Most correla- 
tions were obtained for small scale cylindrical tanks with 
negligible buoyancy effects. In large systems, the buoyancy 
due to thermal stratification generally impedes the motion of 
the colder and denser jet flow and possibly prevents the jet 
flow from reaching the liquid surface. The mixing time 
correlation recently developed by Lin et al. 6 for a liquid 
hydrogen tank included the effect of buoyancy due to the 
thermal stratification. The Richardson number (Ri) was used 
to characterize the buoyancy effect. The study of Lin et al 6 
showed that the mixing time required to destroy the thermal 
stratification increases significantly with the increase in Ri 
when Ri > 0.5. However, the buoyancy effect can be neglected 
if Ri <0.5. The Richardson number was defined based on the 
jet submergence depth (Hj) and the temperature difference 
between the interface and the jet (T s -Tp at the begining of 
mixing process. 

Several experiments 7 ' 11 have been conducted to study 
the steady state vapor condensation rate at the liquid-vapor 
interface in a tank with jet-induced mixing. For high jet 
submergence, Hj/D s >3, Sonin et al. 10 developed a vapor 
condensation rate correlation for a steam-water system with 
negligible buoyancy effect. Brown et al. 1 1 extended this work 
toincludetheeffectofthelowjetsubmergence,0.5 < H:/D s < 3. 
The condensation correlation developed in Refs. 10 and 1 1 
was expressed in terms of liquid subcooling (T s -T b ). The 
correlation was modified by Lin et al. 6 and was expressed as 
a function of jet subcooling (T s -Tj). The modified steady state 
condensation rate correlation was then used to predict the 
transient condensation rate in the fluid mixing of a hydrogen 
tank using the instantaneous value of jet subcooling. The 
limited liquid hydrogen data of Ref. 6 is fairly described by 
this modified condensation rate correlation. 

The experimental work reported in Ref. 6 was performed 
with a small (0. 14 m 3 ) cylindrical tank with relatively low wall 
heat fluxes (3 to 30 W/m 2 ). Therefore, the effect of free 
convection boundary layer (along the tank wall) on the jet- 
induced mixing process was negligible. However, for a typical 
cryogenic storage system, the tank size may be large enough 
such that the flow rate of the warm fluid in the free convection 
boundary layer may be significant. This warm flow moves 
along the tank wall towards the interface and turns inward to 
oppose the subcooled jet. Consequently, evaporation may 
occur at the outer part of the interface (due to free convection 
flow) with condensation occurring at the inner part of the 
interface (due to jet flow). Therefore, the average vapor 
condensation rate at the interface will be reduced and the 
mixing time will be increased. The idealized flow pattern and 
phenomena are shown in Fig. 1. 

This paper presents the experimental results in a large 
storage tank where the liquid hydrogen is mixed by an axial 
turbulent jet directed normally towards the liquid-vapor 
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interface. The effect of jet volume flow rate and liquid fill 
level on the mixing time and the average vapor condensation 
rate at the interface is investigated. The influence of free 
convection flow due to wall heat leak on the mixing process 
and the resulting tank pressure change is also examined. The 
results are compared with those obtained for a small scale 
tank. 6 

The Experiment 

The tests were conducted at the NASA Lewis Research 
Center’s K-Site Facility located at Plum Brook Station in 
Sandusky, Ohio. The test facility had a 7.6 m diameter 
spherical vacuum chamber enclosing a 4.0 m diameter cylin- 
drical shroud which in turn enclosed the LH 2 test tank. The test 
tank was suspended by fiberglass composite struts. The shroud 
was maintained at a constant temperature of 350 K by electri- 
cal resistance heating to obtain a constant heat input to the test 
tank. The pressure in the vacuum chamber during the test 
series was on the order of 10 -4 to 10' 3 Pa. 

The LFL, test tank was insulated with 2 blankets of multi- 
layer insulation, each blanket consisted of 17 layers of double 
aluminumized Mylar separated by silk netting. The tank was 
constructed of 2219 aluminum and was approximately an 
ellipsoidal volume of revolution. The major-to-minor axis 
ratio was 1 .2 with the major diameter of 2.2 m. The test tank 
had a volume of 4.89 m 3 , an internal surface area of 14 m 2 > and 
a mass of 149 kg. Most of the wall was 2.08 mm thick except 
for the thick bolted flange and lid at the top, thickened lands 
for support lugs, and a thickened equatorial region. 

The experimental setup was basically the same as had 
been used for the self-pressurization tests in Refs. 1 and 12. 
However, two systems were added in the present study: a 
pressurization system and a mixing system. The pressuriza- 
tion system was designed such that the pressurant (GH 2 ) was 
supplied from high pressure storage bottles at ambient tem- 
perature. A diffuser was installed near the top of the test article 
and the gas pressurant was fed into the vapor region directly. 
The pressurant flow rate was calculated from pressure drop 
measurements across a square edged orifice placed in the 
pressurant line. The orifice was instrumented with high and 
low range differential pressure transducers as well as an 
upstream pressure transducer and temperature sensor. A mixer 
unit was installed inside the test tank to generate axial jet- 
induced mixing. The mixer unit was the one designed for the 
Shuttle Centaur hydrogen tank and was the same one used in 
Ref. 6. Figure 1 shows a schematic that depicts the mixer as 
installed in the elliptical test tank. The jet nozzle outlet had an 
inner diameter of 0.022 1 m and was located in the center about 
0.51 m from the tank bottom. 

A capacitance probe was installed to measure the liquid 
fill level with an accuracy of ±0.019 m. Instrumentation rakes 
were used for the mounting of silicon diodes temperature 
sensors for temperature measurement. Locations of various 


temperature sensors in the test tank are shown in Fig. 2. The 
fluid temperatures and external tank wall temperature meas- 
urements had an accuracy of ±0.1 and ±0.6 K, respectively. 
Tank pressure was measured by pressure transducers in direct 
communication with the tank ullage with an accuracy of 
±0.7 Pa. The pressurant flow rate measurement had an esti- 
mated accuracy of ±0. 1 8 and ±0.40 kg/hr for the low and high 
range differential pressure transducers, respectively . The mixer 
flow rate was measured by a turbine flowmeter with an 
estimated accuracy of ±0.02 m 3 /hr. 

The heat leak through the tank wall was determined from 
a boil-off test. The procedure was the same as that of Refs. 1 
and 12. With the shroud temperature maintained at 350 K the 
boil-off rate was about 5.6 SCMH (standard cubic meter per 
hour) which yielded an average wall heat flux of 4.2 W/m 2 . 
In each test, the tank liquid was initially saturated at about 
1 17 kPa. The tank was then pressurized either by using the 
pressurization system to introduce GH 2 into the tank or by 
self-pressurization due to the ambient heat leak through the 
tank wall. If the tank was self-pressurized, the bulk liquid 
would be subcooled and thermally stratified. If the tank was 
pressurized by the introduction of GH 2 , thermal stratification 
would be limited to a very thin thermal layer underneath the 
liquid surface. The mixer was turned on when the tank 
pressure reached a specified level. The test was terminated 
when the tank pressure reached a steady value or if the planned 
testing time was reached. 

Results and Discussion 

Three series of experiments consisting of ten tests were 
conducted. The tests were conducted at about 85 and 50 per- 
cent liquid fill levels (in volume) with the jet flow rate ranging 
from 0.91 to 5.63 m 3 /hr. The initial condition of each mixing 
test was established by raising the tank pressure by self- 
pressurization or by the addition of gaseous hydrogen such 
that the bulk liquid was subcooled. Table I shows the experi- 
mental conditions for all the tests. 

Within each test series, the initial conditions such as 
liquid fill and liquid subcooling were very similar and the only 
controlling variable was the jet volume flow rate. Therefore, 
within each test series, the effect of free convection flow 
(along the tank wall) on the mixing of tank fluid induced by 
various jet flow rates can be investigated. Figures 3(a) to (c) 
show the tank pressure decay as a function of time for various 
jet volume flow rates for all the test series. The pressure decay 
rate is faster for a higher jet flow rate as a result of the enhanced 
vapor condensation at the liquid- vapor interface. However, in 
test run 438, which has the lowest jet flow rate, the tank 
pressure increases with time during the mixing. This indicates 
that the net mass transfer at the liquid-vapor interface is 
evaporation due to the free convection flow along the tank 
wall. Also, the jet momentum may not be sufficient to over- 
come the opposing effect of buoyancy due to the thermal 
stratification. 


3 



A solution for a turbulent free convection boundary layer 
with constant wall heat flux and uniform bulk temperature has 
been obtained by Vliet et al. 13 using the velocity and tempera- 
ture profiles employed by Eckert and Jackson 14 in their 
analyses of turbulent free convection on a vertical plate at 
constant wall temperature. Using Vliet et al.’s results, the 
volume flow rate in the boundary layer for a cylindrical tank 
can be written as 

Q bI (x) = 0.292rtDv 1/7 (gpq w /pC p ) 2/7 x 8/7 (1) 

where D is the tank diameter and x is the distance measured 
from the tank bottom. It is obvious that Q bl increases signifi- 
cantly with the tank size. Equation (1) is used to estimate Q bl 
for the nearly ellipsoidal tank of the present study. With the 
modified Grashof number defined as Gr* = g(Jq w Hj 4 / k v 2 , the 
volume flow rate in the boundary layer at the interface location 
is estimated to be 

Q bl * = 0.292 n D s v(Gr*/Pr) 2/7 (2) 

where D s is the diameter of the liquid- vapor interface. Equa- 
tion (2) is expected to over-predict the the volume flow rate 
in the natural convection boundary layer for an ellipsoidal 
tank. This is because of the effects of the tank wall curvature 
and the heat leak through the tank bottom. Due to buoyancy, 
the heating of the tank bottom generally creates a mixing 
region near the tank bottom and delays the formation (shifted 
upwards) of the thermal boundary layer at the sidewall due to 
sidewall heating. 15 

The analysis for a turbulent free circular jet is used to 
estimate the jet volume flow rate at the liquid- vapor interface 
location. According to Schlichting 16 , the volume flow rate 
(Qc), which increases with distance from the orifice due to the 
entrainment of the flow from the surroundings, can be calcu- 
lated by 

Q c (x) = 0.456 Qj ^ (3) 

where Qj is the jet volume flow rate at the jet nozzle outlet, x 
is the distance from the nozzle outlet, and d is the diameter of 
nozzle. The total jet volume flow rate at the liquid-vapor 
interface location is estimated to be 

H i 

Q * = 0.456 Qj— ^ (4) 

c J d 

The values of Q b |* and Qg* are listed in Table n. Fluid 
properties at 1 10 kPa saturated condition were used to calcu- 
late Q bl * 

The jet induced flow brings colder fluid to the central 

part and the boundary flow Q w * brings warmer fluid to the 


outer part of the surface. The net mass transfer at the interface 
results from condensation in the central part and evaporation 
in the outer part of the interface due to the flow rates and 
Q bl * respectively. Also, buoyancy due to thermal stratifica- 
tion may impede the upward motion of the colder and denser 
jet flow and may significantly reduce the condensation rate at 
the liquid- vapor interface. The findings of Refs. 6, 10, and 1 1 
suggest that buoyancy effects may be quantified in terms of 
the Richardson number, Ri. Based on the definition of Ri used 
in Ref. 6, the buoyancy effects due to thermal stratification of 
the bulk liquid are insignificant if Ri < 0.5. Comparison of 
Q bl * and in Table II and the pressure change rate during 
the mixing as shown in Figs. 3(a) to (c) shows that the effect 
of free convection flow along the tank wall on the jet mixing 
process decreases as the jet volume flow rate (Qj) and thus Q c * 
increases. In the test run 438, the jet flow rate Qj is low such 
thatQ c * and Q bl * are comparable. The pressure increases with 
time during the mixing as shown in Fig. 3(a) which suggests 
that the net mass transfer at the liquid-vapor interface is 
evaporation. The condensation rate in the test run 438 may 
also have been significantly reduced due to relatively high 
value of Richardson number (Ri = 0.52, in Table HI). In the test 
run 459, the tank pressure during mixing decreases even 
though Q bI * is greater than Q c *. This is, in part, due to the fact 
that the Richardson number for this test is low (Ri = 0.2, in 
Table HI). If the jet flow rate is high enough such that Q c * is 
much greater than Q^* then there may be condensation 
occuring over most of the interface. As shown in Table n, the 
test run 471 represents this case. 

As shown in Table I, test run 385 has the same jet volume 
flow rate as test run 438. However, unlike test run 43 8, the tank 
pressure in test run 385 decreases (Fig. 3(c)) during the mixing 
process. This is because, for tank pressurization using gaseous 
hydrogen, the pressurization time required to reach the target 
tank pressure is too short for the free convection flow to be 
fully developed and also the liquid temperature is less strati- 
fied. Therefore, the effect of wall heat leak in the mixing 
process is greater in test series A than in test series C. 

The ratio of the liquid-vapor interface diameter to the jet 
nozzle diameter, D g /d, ranged from 85 to 100 in all the tests. 
Thus, as the jet moved upward and reached the liquid surface, 
it had to travel a long distance in the radial direction before 
moving back towards the tank bottom. This made it difficult 
for the axial jet to induce complete mixing of the tank fluid. 
Therefore, unlike the experiments of Ref. 6 in which D $ /d was 
about 25, the thermal equilibrium state was not reached in 
these tests. Each test was terminated when the tank pressure 
reached a steady value or if the planned testing time was 
reached. Table in shows the mixing characteristics for all the 
tests. There are only five test runs (out of ten) in which the tank 
pressure reached a steady value when the test was terminated. 
In Ref. 6, the mixing time was defined as the time required for 
the tank to reach the equilibrium pressure. Thus, the mixing 
time measured for these five test runs can be used to compare 
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the predictions of mixing time correlation developed in Ref. 6 
for a small scale cylindrical tank: 

t u d f 3.3 Ri<0.5 

t * = m j = { (5) 

m D s 2 [2.6 + 2.1 Ri Ri > 0.5 

From Table III, it is seen that the mixing time for test runs 
in which the tank pressure reached a steady value is under- 
predicted by Eq. (5). This may be due to the combined effects 
of natural convection flow, tank size, and tank geometry. The 
effect of natural convection flow can be investigated from the 
present data. Since in each test series, the modified Grashof 
number (Gr*) is nearly constant, the effect of natural convec- 
tion flow on the central jet flow will decrease as the jet 
Reynolds numbers (jet volume flow rates) increases. There- 
fore, as shown in Table HI, the comparison improves as the jet 
Reynolds number (Rep and then the ratio Q c */Q bl * increase. 
Also, the comparison seems to be slightly better for test 
series C because the natural convection flow is not fully 
developed in these tests. 

Vapor condensation at the interface is the key mechanism 
to reduce the tank pressure. Sonin et al. 10 obtained an empiri- 
cal correlation for the rate at which a pure saturated vapor 
condenses on the free surface of a turbulent liquid in a steam- 
water system, under conditions where buoyancy effects are 
insignificant. The vapor was saturated and the vapor conden- 
sation was controlled by the interface turbulence at the liquid 
side. Their vapor condensation rate correlation is expressed as 

st c = 0.0198 Pr -0 ' 33 (1 - 0.5 Ja) (6) 

where St c , defined as St c = m c h fg / [A s p C p v s (T $ -T b )], is the 
average condensation Stanton number , Ja is the Jakob num- 
ber, and m c is the average vapor mass condensation rate at the 
interface. Equation (6) is strictly valid for Ja « 1. The liquid 
subcooling (T $ -T b ) and average turbulence intensity at the 
interface (v g ) are the controlling parameters for the vapor 
condensation rate. 

The average turbulent r.m.s. velocity at the interface, v g , 
for low jet submergence can be determined by 

V S =^ 10.04 p 2 -(7.14 p 2 -3.06 
D s L 

for 0.5 < Hj /D s <3 (7) 

where pj = 0.34 and p 9 = 0.24 from Ref. 17 and Pj = 0.33 and 
P 2 = 0.23 from Ref. 1 8 . For Ja « 1 , the steady state correlation 
of Eq. (6) has been modified by Lin et al. 6 for the application 
of transient cases: 


m c h f g 

S ' c = 

x [l + (D s / d) 2 (0.0198 Pr -0 - 33 ) (v $ / Uj )] 

= 0.0198 Pr -0 ' 33 (1 - 0.5 Ja) (8) 

where the liquid subcooling (T s -T b ) has been replaced by jet 
subcooling (T s -Tp which is usually measured during the test. 
The jet subcooling as a function of time is shown in Figs. 4(a) 
to (c). According to the range of jet subcooling shown in 
Figs. 4(a) to (c), the Jakob number is less than 0.05 for all the 
tests. 

Equation (8) is used to predict the change of vapor mass 
for the tests in test series A. The reason to use test series A for 
comparison is that Hj/D s is about 0.5 and the vapor superheat 
for these tests is relatively mild. The vapor is superheated with 
the temperature difference (T dome -T s ) ranging from 40 to 
100 K for test series A. These conditions are relatively close 
to the applicable range of the correlation. The vapor mass is 
obtained from a curve-fitted one-dimensional temperature 
distribution. This one-dimensional calculation is a good 
approximation because the data indicates that the radial varia- 
tion of the vapor temperature is small. 

Figures 5(a) to (c) show the estimated vapor mass as a 
function of time. The vapor mass predicted by using Eq. (8) is 
also included in the figures for comparison. According to the 
previous discussion, Eq. (8) should over-predict the average 
condensation rate at the liquid-vapor interface because of the 
effect of natural convection flow. Figure 5(a) indeed shows 
this fact. As the jet volume flow rate increases, the comparison 
gets better. Surprisingly, Figs. 5(b) and (c) show that Eq. (8) 
under-predicts the average condensation rate. This can be 
explained by the fact that, in Figs. 5(b) (for test run 434) and 
(c) (for test run 471), the jet volume flow rate is sufficiently 
high that surface breakup occurs. The agitation of the liquid- 
vapor interface induces significant condensation and Eq. (8), 
which was developed based on the assumption of a flat 
interface, does not apply. In Ref. 7, an equation was provided 
to estimate the condition for the onset of surface breakup: 

Rej >0.46(Apga) 1/4 Hj/(p 1/2 v) (9) 

where Ap is the density difference between the liquid and the 
vapor and o is the surface tension of the liquid. Equation (9) 
gives the prediction of minimum jet Reynolds number (Rep 
for the jet to penetrate the liquid- vapor interface. Using this 
criteria, the jet Reynolds number for the onset of the surface 
breakup for the tests at 85 percent fill level is about 230 000. 
Comparing with the jet Renolds numbers listed in Table HI, 
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the jet Reynolds numbers for both test runs 434 and 471 are 
greater than 230 000. This comparison supports our argument 
associated with Figs. 5(b) and (c). 

Conclusion 

Experiments have been conducted to investigate the 
effect of fluid mixing on the depressurization of a large scale 
liquid hydrogen storage tank. Under the experimental condi- 
tions in the present study, the following conclusions can be 
drawn: 

1. For a large cryogenic fluid storage tank subjected to 
very low wall heat flux of 4.2 W/m 2 , the effect of natural 
convection flow on the jet induced mixing process may not be 
neglected. The boundary layer flow of warm fluid (along the 
tank wall) due to natural convection decreases the average 
condensation rate at the liquid-vapor interface and increases 
the mixing time. 

2. The dimensionless mixing time for the tests under 
comparable conditions (Table El) is under-predicted by the 
correlation (Eq. (5)) developed based on small scale tank 
results. 

3. The thermal equilibrium state is very difficult to 
achieve for fluid in a large scale mixing tank, i.e., in a tank with 
large ratio of liquid surface diameter to jet nozzle diameter. 
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TABLE L— EXPERIMENTAL CONDITIONS FOR MIXING TESTS 


Test 

run 

no. 

Pressurization 

source 

Initial Liquid Fill 

(m 3 Air) 

D s 

(m) 

(kPa) 

p f 

(kPa) 

(Vty 

(K) 

% 

V,(m 3 ) 

Test Series A 

438 

S 

86.9 

4.25 

0.91 

1.87 

187.0 

HI 

1.7 

436 

S 

85.3 

4.17 

1.82 

1.91 

187.0 

- n 

1.7 

434 

S 

86.3 

4.22 

3.47 

1.88 

186.1 


1.7 

471 

S 

81.6 

3.99 

5.63 

1.99 

186.3 


1.9 

Test Series B 

459 

S 

48.3 

2.36 

0.91 

2.22 

186.3 

158.3 

1.6 

457 

S 

49.1 

2.40 

1.82 

2.22 

186.5 

149.5 

1.6 

449 

S 

49.1 

2.40 

3.41 

2.22 

186.1 

145.8 

1.6 

Test Series C 



85.3 

B9 

0.93 

19 

Ml 

158.7 

1.7 



84.3 

W 

1.82 

19 

ESI 

127.8 

1.7 

389 


84.9 

HI 

3.40 

19 

19 

121.8 

1.7 


S: Self-pressurization G: Gaseous hydrogen pressurant 


TABLE 2 — COMPARISON OF VOLUME FLOW RATE AT THE 
INTERFACE LOCATION DUE TO THE CENTRAL JET AND 
THE NATURAL CONVECTION FLOW 



Jet submergence 

Qc* 

Liquid Height 

Interface 

Gr* 

Qbi* 


depth, Hj(m) 

(m 3 /hr) 

Hj(m) 

Dia., D s (m) 


(m 3 /hr) 

Test Series A a 

438 

Ha 

17.8 


1.87 

9.2X10 14 

20.0 

436 

1 

34.5 


1.91 

8.4X10 14 

20.0 

434 

0.94 

67.3 


1.88 

8.9X10 14 

20.0 

471 

0.86 

99.9 


1.99 

7.1X10 14 

19.8 

Test Series B a 

459 

0.41 

7.7 

0.92 

2.22 


14.0 

457 


15.8 

0.93 

2.22 


14.2 

449 

9mm. 

29.6 

0.93 

2.22 


14.2 

Test Series C b 

— 

385 

0.92 

17.7 


um 


m 

387 

0.90 

33.8 


MEM 



389 

0.91 

63.8 


MEM 




a Self-pressurization. 

b Gaseous hydrogen pressurant. 

c The natural convection flow is not fully developed. 
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TABLE 3 — THE MIXING CHARACTERISTICS OF THE TESTS 


Test 
Run No. 

Re ; 

Ri 

H/D s 

u j 

(m/s) 

t * a 
m 


Qc*'Q b ,* 

Test Series A 

438 

80 000 

0.52 

0.51 

0.66 

— 

— 

0.89 

436 

160 000 

0.13 

0.48 

1.32 

13.1 c 

3.3 

1.73 

434 

304 000 

0.036 

0.50 

2.51 

7.49° 

3.3 

3.37 

471 

495 000 

0.014 

0.43 

4.08 

9.11 c 

3.3 

5.05 

Test Series B 

459 

80 000 

0.20 

0.18 

0.66 

15.1 

3.3 

0.55 

457 

161000 

0.049 

0.19 

1.33 

12.7 

3.3 

1.11 

449 

299 000 

0.014 

0.19 

2.47 

16.8° 

3.3 

2.08 

Test Series C 

385 

82 000 

0.49 

0.48 

0.68 

2.3 

3.3 

— 

387 

161 000 

0.13 

0.47 

1.33 

4.4 

3.3 

— 

389 

296 000 

0.038 

0.47 

2.44 

8^ 

3.3 

— 


a Dimensionless urae based on the measurement- 

^Dimensionless time based on the prediction of Eq. (5). 

c Test was terminated when the tank pressure reached a steady value. 



Qw 

Figure 1 . — The flow pattern and phenomena of a mixing tank. 


• Wall temperature sensor 
■ Fluid temperature sensor - operating 
□ Fluid temperature sensor - failed 
| Interface rake - 0.64 cm spacing 
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Figure 3. — Pressure decay as a function of time, (a) For test series Figure 4. — Jet subcooling as a function of time, (a) For test series A. 

A. (b) For test series B. (c) For test series C. (b) For test series B. (c) For test series C. 
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Figure 5. — The change of vapor mass as a function of time, (a) For 
test run 436; Qj = 1 .82 m 3 /hr. (b) For test run 434; Qj = 3.47 m 3 /hr. 
(c) For test run 471 ; Qj * 5.63 m 3 /hr. 


10 






REPORT DOCUMENTATION PAGE 


Form Approved 
OMBNo. 0704-0188 


Public reporting burden for this collection of information is estimated to average t hour per response, incfurSng the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302. and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY ( Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 

June 1994 


3, REPORT TYPE AND DATES COVERED 

Technical Memorandum 


Experimental Investigation of Jet-Induced Mixing of a Large Liquid Hydrogen 
Storage Tank 


6. AUTHOR(S) 


C. S. Lin, M. M. Hasan, and N. T. Van Dresar 


5. FUNDING NUMBERS 


WU-506-42-73 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-8926 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONfTORING 
AGENCY REPORT NUMBER 

NASA TM- 106629 
AIAA-94— 2079 


11. SUPPLEMENTARY NOTES 

Prepared for the 6th Joint Thermophysics and Heat Transfer Conference cosponsored by the American Institute of Aeronautics and Astronautics and the 
American Society of Mechanical Engineers, Colorado Springs, Colorado, June 20-23, 1994. Chin-Shun Lin, Analex Corporation, 3001 Aerospace 
Parkway, Brook Park, Ohio 44142 (work funded by NAS3-25776); Mohammad M. Hasan and Neil T. Van Dresar, NASA Lewis Research Center. 
Responsible person, Mohammad M. Hasan, organization code 5340, (216) 433-7494. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Category 34 


13. ABSTRACT (Maximum 200 words) 

Experiments have been conducted to investigate the effect of fluid mixing on the depressurization of a large liquid 
hydrogen storage tank. The test tank is approximately ellipsoidal, having a volume of 4.89 m 3 and an average wall heat 
flux of 4.2 W/m 2 due to external heat input. A mixer unit was installed near the bottom of the tank to generate an upward 
directed axial jet flow normal to the liquid- vapor interface. Mixing tests were initiated after achieving thermally strati- 
fied conditions in the tank either by the introduction of hydrogen gas into the tank or by self-pressurization due to 
ambient heat leak through the tank wall. The subcooled liquid jet directed towards the liquid-vapor interface by the 
mixer induced vapor condensation and caused a reduction in tank pressure. Tests were conducted at two jet submergence 
depths for jet Reynolds numbers from 80 000 to 495 000 and Richardson numbers from 0.014 to 0.52. Results show that 
the rate of tank pressure change is controlled by the competing effects of subcooled jet flow and the free convection 
boundary layer flow due to external tank wall heating. It is shown that existing correlations for mixing time and vapor 
condensation rate based on small scale tanks may not be applicable to large scale liquid hydrogen systems. 


14. SUBJECT TERMS 


Liquid hydrogen storage tank; Jet mixing; Pressure decay; Vapor condensation 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


15. NUMBER OF PAGES 

12 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 























